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1.  SUMMARY 

The  results  of  the  second  year  of  this  prograa  include:  (1)  A  complete 
re-evaluation  of  Raasan  data  previously  taken  in  an  Res8500  turbulent 
nonpreaixed  syngas  flaae,  shoving  significant  corrections  due  to  high- 
teaperature  effects,  (2)  a  detailed  coaparison  vith  Monte-Carlo/partial- 
eqni librium  calculations,  bolstering  the  concept  of  partially-equilibrated 
radical  pools  in  turbulent  coabustion,  (3)  the  developaent  of  CO/I^/l^  mix- 
tures  vith  lov  extinction  liaits  (approaching  hydrocarbons)  vith  Rayleigh- 
scattering  cross-sections  aatched  for  density  and  vith  viability  for  Raaan 
aeasureaents,  (4)  a  broad  range  of  iaproveaents  to  ainiaixe  interference  due 
to  cheai-luaine sconce  in  Raaan  aeasureaents,  (5)  the  design  and  construction 
of  a  co- annular  preaixed  pilot  for  the  nonpreaixed  burner,  (6)  detailed  Raaan 
aeasureaents  at  Re~5000  and  Re«15,000  in  highly-strained  regions  as  veil  as 
further  dovnstreaa  in  the  nev  pilot-stabilised  flaae,  (7)  coaparison  vith  a 
strained  flsaelet  theory  for  localised  extinction  and,  (8)  the  developaent  of 
a  theory  vhich  recognises  that  a  vide  range  of  Daakohler  nuabers  exists  and  so 
treats  only  fast  reactions  in  the  context  of  flaaelets.  The  latter  theory 
dravs  froa  earlier  vork  on  partial-equil ibr iua.  Taken  together,  these  results 
are  a  step  towards  understanding  the  nature  of  localised  turbulent  estinction, 
a  precursor  of  blowout  in  practical  devices  under  soae  conditions.  In  parti¬ 
cular,  the  popular  notion  of  thin  flaaelets  as  described  by  high-energy  asyap- 
totics  is  in  serious  difficulty,  because  of  the  vide  spectrua  of  Daakohler 
nuabers  associated  vith  the  aulti-atep  nature  of  cheaical  kinetics  of  coabus- 
tion. 
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2.  RESEARCH  OBJECTIVES 

The  research  objectives  of  this  program  are  to: 

A.  Iaprove  and  expand  the  database  for  turbulent  carbon 
aonoxide/hydrogen/nitrogen  jet  flaaes  utilizing  pulsed  Ran an  spectros¬ 
copic  aeasureaents. 

B.  Use  experiaental  results  to  evaluate  assuaptions  Bade  in  foraulating 
prediction  aethods  based  on  probability  denaity  functions  and  on  laainar 
flaaelets. 

C.  Modify  the  experiaental  facility  to  oonduct  turbulent  jet  coabustion 
tests  at  Reynolds  nuabers  greater  than  9,000  without  blow-off  at  nozzle. 
Use  PL IF  to  visualize  extent  of  extinction. 

D.  Model  a  carbon  aonoxide/ hydrogen/nitrogen  jet  flaae  at  S,500  Reynolds 
nuaber  using  two-scalar  and  flaaelet  approaches. 

E.  Conduct  experiaents  on  two  carbon  aonoxide/hydrogen/nitrogen  jet  flaaes 
at  Reynolds  nuabers  above  9,000  under  conditions  of  local  extinction. 

F.  Foraulate  alternatives  to  the  partlal-equil ibr iua  cheaical  kinetics 
aodel  previously  used  to  predict  the  cheaistry  of  carbon  aonoxide 

flaaes. 


3.  STATUS 


3.1  Introduction 

Development  of  a  fundamental  understanding  of  turbulence-chemistry 
interactions  is  one  of  the  most  important  and  challenging  problems  in  tur¬ 
bulent  combustion.  Although  the  fast  chemistry  or  "mixed- is-burned"  model  for 
i  action  is  useful  for  understanding  gross  aerodynamics  in  many  combustion 
systems,  there  is  an  increasingly  important  range  of  phenomena  that  cannot  be 
so  accounted  for.  For  example,  finite-rate  radical  pool  chemistry  has  been 
shown  to  be  very  influential  in  thermal  NO^  fclu|  CO  emissiona  from  atmospheric 

and  CO/Hj  flames  [Drake  et  al.  (1987)]  and  to  degrade  combustion  efficiency 
in  supersonic  combustion  systems  [Correa  and  Mani  (1987)].  Localised  effects 
including  extinction  due  to  turbulent  stretching  may  be  largely  dne  to  the 

chemistry  not  being  fast  compared  with  the  turbulent  microscales  [Masri  and 

Bilger  (1984)].  as  has  been  demonstrated  in  laminar  CB^  flames  [Tsuj  i  and 

Yamaoka  (1970)].  High-altitude  limitations  on  aircraft  combustors  are  caused 

in  part  by  slow  kinetics  due  to  the  lorn  densities.  It  is  clear  that  as  the 
limits  of  combustion  equipment  are  expanded,  it  is  more  and  more  important  to 
account  for  finite-rate  chemical  kinetics.  On  the  other  hand,  detailed  chemi¬ 
cal  kinetic  mechanisms  have  been  established  only  for  the  combustion  of  rela¬ 
tively  simple  fuels  such  as  CO/^  and  CH4  [Varnats  (1981)].  The  combus¬ 
tion  of  complex  hydrocarbons  of  practical  significance  is  usually  represented 
by  a  semi-empirical  first  step  to  simple  constituents  such  as  CO  and  Hj,  with 
subsequent  oxidation  of  the  latter  (Varnats  (1981)].  This  may  compromise  the 
understanding  of,  for  example,  "prompt"  NO^  wgiCg  is  formed  via  hydrocarbon 
fragments  [Iverach  et  al.  (1972)] t  it  also  leaves  only  the  relatively  simple 
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fuels  as  candidates  for  fundamental  work  in  tnrbnlent  flames,  where  the  chem¬ 
istry  must  be  known  a  priori. 

The  approach  of  this  joint  ezperimental/computational  program  is  to  study 
turbulence-chemistry  interactions  in  pilot-flame  stabilized,  turbulent,  non- 
premized  jet  ("diffusion")  flames  of  Hj  and  CO/Hj/Nj  fuels,  where  the  tur¬ 
bulent  fuel/air  mizing  and  combustion  chemistry  are  reasonably  well  known. 
(The  turbulence  is  understood  better  than  that  in  the  recirculating  flowfields 
characteristic  of  practical  burners).  Ezper imental ly  our  unique  combustor 
facilities  and  laser  measurement  teohniques  are  prowiding  time-  and  space- 
resolved  pointwiae  meaaurementa  of  probability  density  functions  of  tempera¬ 
ture,  major  apecies  concentrationa,  density  and  miztnre  fraction  in  the  flow. 
Two-dimensional  images  of  the  density  field  using  planar  Kayleigh  scattering 
will  indicate  the  local  stretch  (instantaneous  gradient  of  compositional 
scalar).  Ezperimental  results  on  stretched  laminar  opposed-flow  diffusion 
flames  are  also  awailable  for  comparison.  Similar  images  of  the  OH  field  will 
indicate  the  eztent  of  localised  eztinction.  Computationally,  the  single¬ 
scalar  fast  chemistry  model  with  assumed  shaped  pdf,  joint  pdf  models  for 
scalar  and  dissipation  rate  closed  with  the  laminar  flamelet  model,  [Williams 
(1974)]  and  two-  or  three-scalar  non-equil ibrium  approaches  using  mizture 
fraction  and  reaction  progress  variables  are  being  developed.  Detailed 
ezper iment- model  comparisons  are  being  used  to  assess  models  and  so  to  provide 
a  more  fundamental  understanding  of  turbulence-chemistry  interactions  and 
eztinction. 
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3.2  Turbulent  Flue  Facility:  Summary  of  Previoua  Work 

Laser  diagnostic  techniques  for  planar  OB  fluorescence,  laser  velo- 
ciaetry,  pulsed-Ruan  scattering,  and  laser-saturated  OH  fluorescence  have 
been  used  in  our  laboratory  to  characterize  turbulent  jet-dif fuaion  flaaes 

^  an(l  CO/H2/N2  fuels.  The  characteristics  of  the  unstabilized  coabustor 
and  the  initial  conditions  for  the  flues  have  been  discussed  in  detail 
[Correa  et  al.  (1984),  Lapp  et  al.  (1983)].  The  coabustor  consists  of  a  3.2 
aa  internal-diueter  fuel  tube  centered  axially  in  a  15-ca- square  by  1-a  long 
test  section.  The  fuel  velocity  was  controlled  by  calibrated  critical  flow 
orifices  and  the  velocity  of  the  surrounding  co-flowing  air  streu  by  a 
servo-control  on  the  exhaust  fan.  For  the  H^  flues,  the  fuel  jet  had  an  ini¬ 
tial  average  velocity  of  283  a/s  and  is  expected  to  have  a  turbulent  pipe-flow 
profile  corresponding  to  a  cold-flow  Keynolds  nsaber  of  8300.  The  inlet  air 
velocity  was  12.3  a/s  and  was  flat  to  within  ±2%.  The  axial  pressure  gradient 
(deterained  by  14  wall-pressure  taps  and  by  free-streu  velocity  aeasureaents) 
**•  ”51  P»/a.  For  the  CO/Hj/Nj  flue,  the  initial  aolar  ratio  of  the  fuel  was 
approxiaately  40%CO/30%H2/30%h^  with  a  uall  uount  of  CH^  (-  0.7%)  soaetiaes 
added.  The  initial  fuel  and  air  velocities  were  34.6  a/s  and  2.4  a/s,  respec¬ 
tively,  and  were  chosen  to  natch  the  Reynolds  nuaber  (8300)  and  initial  fuel- 
to-air  velocity  ratio  (22.8)  of  the  ^  fi**,.  The  initial  jet  velocity  pro¬ 
file  was  Measured  1  an  froa  the  nozzle  exit,  and  the  axial  pressure  gradient 
was  -1.0  Pa/a. 

The  experiaental  aeasureaents  on  the  turbulent  H^  flues  are  sunaarized 
elsewhere,  particularly  in  Drake  et  al.  (1986).  The  data  include  qualitative 
aad  quantitative  iaaging  (Schlieren,  shadowgraph,  and  planar  OH  fluorescence ) j 


instantaneous  measurements  of  velocity  (laser  vel ocinetry) ;  siaaltaneoos 
aeasoreaents  of  teaperatore,  density,  individual  major  speciea  concentrations, 
and  aixtnre  fraction  (pul sed-Raaan) j  and  instantaneoos  aeasoreaents  of  OH  con¬ 
centrations  (laser-saturated  fluorescence).  A  collection  of  instantaneous 
aeasureaents  at  a  given  flaae  location  provide  probability  density  functions 
(both  conventional  and  Favre  averaged).  Analysis  of  those  pdf's  have  permit¬ 
ted  determination  of  interaittency  and  conditional  moments  [Pits  and  Drake 
(1986)].  The  first  four  aoaents  of  the  aizture  fraction  pdf  show  reaarkable 
siailarity  to  the  aoaents  aeasuxed  in  turbulent  nonreacting  jet  flows  and 
indicate  the  presence  of  broad  airing  zones  [Drake  et  al.  (1985)].  Coaparison 
of  the  Kaman  data  with  thermodynamic  and  fluid  mechanic  models  quantify  the 
importance  of  preferential  diffusion  [Drake  et  al.  (1982)]  of  and  0f  radi¬ 
cal  super equil ibrium  [Drake  et  al.  (1986)].  Superequil ibriua  was  confined  by 
direct  aeasureaents  of  OH  concentrations  using  laser-saturated  OH  fluorescence 
with  aeasured  mean  OH  concentrations  far  above  adiabatic  equil ibriua  values 
close  to  the  nozzle  and  near  equilibriua  values  far  downstreaa. 

The  ezperiaental  aeasureaents  in  the  turbulent  CO/H^/l^  flaae  are  surnaar- 
ized  in  Correa  et  al .  (1984)  and  Drake  (1985).  The  data  include  velocity, 
teaperature,  density,  aizture  fraction,  and  individual  species  concentrations 
as  in  the  H^  flame.  In  addition,  probe  saapling  aeaaureaents  of  average  value 
of  NO  and  NO^  mere  obtained  for  the  base  fuel  (for  theraal  N0X) ,  with  saall 
aaounts  of  NH^  added  (for  FBN  conversion),  and  with  saall  aaounts  of  CH^  added 
(for  proapt  NO^). 

Extensive  comparisons  with  a  partial-equil ibriua  nodal  for  the  radical 
pool  and  k-s/assuaed  shape  pdf  nodal  for  turbulent  flow  were  reported  by 
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Correa  et  al.  (1984))  a  more  sophisticated  theory  vas  presented  by  Correa 
(1986a).  These  theories  successfully  predicted  the  high  levels  of  supere- 
quilibriua  free  radicals  and  were  extended  to  thenal  NO^  v i*  the  Zeldovich 
aechanisa  in  Drake  et  al .  (1987).  Additionally  Pope  and  Correa  (1986)  used 
partial-equil ibr iua  and  a  Monte-Carlo/ joint  pdf  flo*  siaulation  to  show  that, 
in  the  saae  flaae,  aost  aaj or  species  were  veil  predicted  except  for  CO.  This 
discrepancy  has  been  found,  as  discussed  belov,  to  stea  from  errors  in  the 
reduction  of  the  Raaan  data.  In  fact,  the  partial  equilibriua  model  seeas  to 
be  aore  accurate  than  concluded  by  Pope  and  Correa  (1986).  Correa  et  al. 
(1987)  have  presented  these  argnaents  in  detail. 

In  the  first  year  of  the  current  prograa  [Correa  (1986b)]  difficulties 
associated  with  aeasureaent  of  CO^  by  Raaan  scattering  were  addressed,  the 
laainar  flaaelet  approach  vas  assessed  by  exaaination  of  the  Raaan  data  for 

turbulent  and  CO/Hj/J^  jet  flaaes  in  cofloving  air  at  Reynolds  nuabers  of 
8S00  or  less  and  a  pilot-stabilized  burner  for  high  Reynolds  nuaber  flaaes  vas 
designed.  As  a  consequence  of  these  studies,  the  00^  measurement  in  our  Raaan 
system  is  systeaatical ly  improved)  ve  also  raised  serious  questions  regarding 
the  validity  of  viewing  nonpreaixed  turbulent  flaaes  as  enseables  of  stretched 
laainar  flaaelets.  The  argnaents  will  not  be  reproduced  here. 

Work  performed  in  the  second  year  follows  along  the  saae  lines.  New  data 
and  aodel ing  in  high  Reynolds  nuaber  nonpreaixed  turbulent  flaaes  are  provi¬ 
ding  a  basis  for  evaluating  current  ideas  on  turbulence-cheaistry  interactions 
including  extinction  and  for  formulating  new  approaches. 
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3.3  Reanalysis  of  Prior  Data  and  Comparison  with  Model 

A  significant  improvement  in  the  analysis  of  the  rav  Raman  data  for  the 
Re  =  8500  flame  of  40%  CO,  30%^,  30%  fuel  was  made  by  correcting  for  the 
fraction  of  Raman  signal  in  the  Raman  band  pass  of  the  major  species  (0^,  , 

^2 <  0,  CO)  at  elevated  temperatures.  This  correction  factor  vas  inadver¬ 

tently  set  to  unity  in  prior  analyses  of  the  rav  Raman  signals.  The  signifi¬ 
cance  of  this  correction  to  the  measured  Raman  signal  at  elevated  temperatures 
arises  from  the  fact  that  the  Raman  exit  slits  have  a  finite  width  due  to 
space  limitations  and  noise  considerations)  therefore  significant  portions  of 
the  hot  Raman  bands  at  higher  temperatures  do  not  fall  within  the  finite  exit 
slit  and  are  not  measured.  The  correction  factors  required  to  account  for  the 
limited  Raman  spectral  band  pass  are  a  function  of  the  species  involved,  tem¬ 
perature,  width  of  the  exit  slit  and  the  wavelength  region  covered  by  the  phy¬ 
sical  Raman  band  pass.  They  were  calculated  by  determining  the  fraction  of 
the  predicted  Raman  vibrational  contour  in  the  Raman  bandpass  using  codes 
developed  by  Lapp  et  al.  (1983).  Table  1  lists  typical  values  of  the  correc¬ 
tion  factors  for  the  major  species  involved.  As  can  be  seen,  the  correction 
factors  for  and  HjO  are  significantly  in  excess  of  unity,  whereas  those  for 

*&4  ^2  tr#  !**•  than  unity.  Table  2  compares  typical  mean  values  of 

temperature,  mixture  fraction  and  mole  fractions  of  major  species  obtained  at 
x/da25  in  the  Res8500  syngas  flame  for  two  cases  (a)  all  correction  factors 
set  to  unity  and  (b)  correction  factors  allowed  to  assume  proper  values 

(listed  in  Table  1).  The  mean  temperature  values  are  less  sensitive  to  non¬ 

unity  correction  factors  since  temperature  is  based  upon  the  sor-of-mole- 
fraction  method  which  is  not  significantly  affected  by  the  corrections. 


Table  (1) 


High  temperature  correction  factors  for  vibrational  Raxian  scattering 
osed  to  correct  for  fraction  of  hot  bends  of  various  species  falling 
outside  their  respective  exit  slits. 


Temp. 

( *  K ) 

°2 

CO 

n2 

300 

1.00 

1.00 

1.00 

1.000 

600 

.975 

.994 

1.019 

1.001 

900 

.914 

.966 

1.009 

1.011 

1200 

.848 

.922 

.982 

1.058 

1500 

.799 

.876 

.953 

1.141 

1800 

.769 

.839 

.932 

1.270 

2100 

.756 

.812 

.922 

1.440 

2400 

.757 

.796 

.922 

1.650 

2700 

.768 

.790 

.932 

1.900 

3000 

.787 

.792 

.950 

2.159 

A.  A- 
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Temperature  and  mole  fraction  of  major  species  at  x/d=25  obtained 
without  correcting  for  high  temperature  effects  (correction  factors 

set  to  unity) . 


(30b  1^, 

40%  CO 

.  30% 

Nj )  = 

8500. 

x/d  *  25 

Temperature 

Mole 

Fraction 

Mixture 

Fraction 

Position 

T/2000K 

co2 

°2 

CO 

N2 

®2° 

*2 

Conv. 

Favre 

(r/a) 

.63 

.600 

.084 

.007 

.233 

.489 

.105 

.075 

.443 

.448 

1.89 

.687 

.099 

.012 

.189 

.530 

.114 

.049 

.385 

.39 

3.15 

.722 

.110 

.034 

.127 

.587 

.109 

.024 

.297 

.293 

4.39 

.647 

.101 

.084 

.054 

.655 

.090 

.008 

.199 

.176 

5.66 

.478 

.061 

,139 

.021 

.706 

.062 

.002 

.117 

.081 

6.92 

.333 

.035 

.175 

.010 

.733 

.037 

.001 

.059 

.033 

8.18 

.229 

.017 

.194 

.004 

.7  57 

.019 

.001 

.023 

.011 

9.44 

.159 

.005 

.208 

.001 

.770 

.007 

0.000 

.004 

.001 

10.69 

.146 

.004 

.209 

.001 

.772 

.004 

0.000 

.001 

.000 

'  4*  . 

w. 

X'-s 


.OvO’y&A 


i 


I 


Table  (2b) 

Temperature  and  sole  fraction  of  major 
for  high  temperature  vibrational  bands 


species  at  x/d=25  corrected 
using  the  correction  factors 


1 isted  in  Tabl e  1 . 


Teaiperatur  e 

Mole  Fraction 

Mixture 

Fraction 

Position 

(r/a) 

T/2000K 

C02 

°2 

CO 

n2 

1^0 

*2 

Conv . 

Favre 

0.63 

.638 

.080 

.005 

.193 

.505 

.122 

.090 

.541 

.541 

1.89 

.739 

.093 

.007 

.144 

.545 

.141 

.062 

.504 

.506 

3.15 

.785 

.104 

.025 

.085 

.605 

.141 

.032 

.408 

.396 

4.39 

.699 

.096 

.072 

.025 

.676 

.113 

.010 

.267 

.225 

5.66 

.504 

.059 

.129 

.004 

.725 

.072 

.003 

.145 

.093 

6.92 

.343 

.034 

.169 

0.000 

.746 

.040 

.001 

.067 

.036 

8.18 

.231 

.016 

.191 

0.000 

.763 

.019 

.001 

.024 

.011 

9.44 

.159 

.005 

.207 

0.000 

.771 

.007 

0.000 

.004 

.001 

10.69 

.146 

.004 

.209 

0.000 

.773 

.004 

0.000 

.001 

.000 
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The  result*  of  the  partial-equil  ibr ium/Monte-Carl o  nodel  [Pope  and  Correa 
(1986)]  were  compared  with  data  on  the  noo-premixed  40%  CO,  30%  H^,  30b 
"syngas"  jet  flaaie  in  co-flowing  air.  The  Raman  data  on  major  species  density 
and  temperature  are  those  of  Lapp  et  al.  (1983)  corrected  here  for  high  tem¬ 
perature  effects  as  described  above.  The  Raman  data  are  compared  at  three 
axial  locations,  (x/d*10,  25  and  50) j  x/d=50  corresponds  approximately  to  the 
visible  length  of  the  flame. 

Radial  profiles  of  the  Favre-averaged  mean  mixture-fraction  (£)  agree 
quite  well  at  x/d=10  (Fig.  1)  except  that  the  jet  is  predicted  to  have  decayed 
too  slowly  on  the  centerline,  and  ?  is  slightly  overspread  at  the  boundary. 
These  differences  are  exaggerated  in  the  mean  density  normal  ixed  by  that  of 
air  (Fig.  1).  The  central  region  of  the  jet  is  hotter  (lower  density)  than 
indicated  by  the  data.  Temperature  profiles  (Fig.  2)  further  confirm  this 
discrepancy  and  also  iudicate  a  difference  of  about  1001  in  the  off-axis 
peaks.  Nitrogen  concentrations  plotted  in  Fig.  2  show  better  agreement.  It 
should  be  recalled  that  the  mixture  fraction  data  are  based  on  the  measured 

and  H^o  concentrations.  More  will  be  said  on  this  below.  Some  of  the  other 
mean  major  species  concentrations  are  shown  in  Fig.  3.  Compared  with  the 
data,  CO  is  underpredicted  by  12%  at  worst,  i,  overpredicted  by  about  10%, 
and  C02  overpredicted,  which  is  consistent  with  a  carbon  balance.  The  jet 
spread  is  again  slightly  overpredicted.  The  discrepancies  in  the  core  seem 
consistent  with  overpredicting  the  amount  of  CO  burned  and  therefore  the  tem¬ 
perature.  The  oxygen  profiles  also  appear  to  agree  reasonably  well. 

Mean  mixture  fraction  and  density  at  x/d-25  show  that  the  jet  spread  is 
predicted  quite  well  but  the  centerline  mixture  fraction  is  about  10%  too  low 


(Fig.  4).  Correspondingly,  since  the  gas  is  predicted  to  be  closer  to 
stoichioaetr ic,  the  teaperature  is  higher  by  aboot  150*K  than  aeasured  in  the 
central  region  (Fig.  5).  The  spread  rate  appears  quite  reasonable.  Once 
again,  the  nitrogen  profiles  (Fig.  5)  agree  veil  except  in  the  central  region, 
i.e. ,  r/a<4. 

Mean  major  species  concentrations  at  x/d-25  are  shown  in  Fig.  6.  The 
agreeaent  is  auch  better  than  presented  in  Pope  and  Correa  (1986)  because  of 
the  changes  in  the  Kaaan  data.  Here  CO  is  underpredicted  by  at  vorst  10b  and 

^2*  ^2  an<*  *r#  predicted  very  veil.  In  fact,  the  error  in  CO  is  rather 

aore  suggested  of  an  over-dif fnaive  calculation,  although  of  course  that  is 
not  inherent  vithin  the  Monte-Carlo/ Joint  pdf  aodel,  than  a  systematic  error 
which  aight  be  attributed  to  a  breakdown  of  partlal-equil  ibr  iua  in  the 
radical-pool . 

Mean  aixture  fraction  profiles  at  x/d*30  show  significant  disagreeaent 
whereas  the  aean  density  (Fig.  7)  agrees  well  with  the  data.  Mean  teaperature 
and  concentrations  also  agree  fairly  well  as  seen  in  Fig.  8.  The  good 
*Sr*9a*Ilt  ^  fraction  is  at  first  glance  surprising  given  the  relative 

lack  of  agreeaent  in  £.  is  essentially  an  inert  species  -  only  trace  quan¬ 
tities  are  converted  to  NO^  -  »nd  therefore  should  follow  a  conserved  scalar 
distribution. 

This  apparent  contradiction  can  be  resolved.  The  aeasureaent  of  the  aole 

fraction  of  i»  independent  of  the  systeaatic  errors  involved  in  the 

aeasureaent  of  other  epeciee,  whereas  the  aixture  fraction  is  a  derived  quan¬ 
tity  which  involves  errors  in  aeasureaents  of  all  aajor  species.  The  lnstsn- 


16 


taneous  value  of  fraction,  can  e.g.  be  defined  as 


-  *S2l,'*.u 


Hence,  the  systematic  errors  in  (  are  dne  to  (a)  errors  in  measurement  of  all 
major  species  and,  (b)  neglecting  all  radicals  and  other  minor  species  since 
they  are  not  measured  by  Kaaan  scattering.  These  errors  affect  the  "measured" 
molecular  weight  ^  used  in  Bq.  (1).  In  fact,  a  change  in  mean  molecular 

weight  of  the  gas  mixture  from  27.86  to  28.8,  for  the  same  measured  mole  frac¬ 
tion  of  Nitrogen  (1^  >0.711)  results  in  the  mixture  fraction  changing  from 

0.15  to  0.22.  Thus,  "measurements"  of  the  mixture- fraction  are  very  sensitive 
to  the  molecular  weight  of  the  mixture  especially  at  low  values  of  f  (as  at 
x/d-50) .  It  should  be  noted  that  all  the  data  for  ?  are  based  on  Hydrogen 
element,  but  the  same  argument  holds  since,  as  Table  3  shows,  the  value  of  ? 
based  on  three  elements  (C,H  and  N)  agree  well.  This  result  is  expected  since 
the  effects  of  differential  diffusion,  observed  in  flames  at  lover  Reynolds 
numbers  [Drake  et  al .  (1982)]  should  be  negligible  in  this  (high  Reynolds  num¬ 
ber)  flame.  Systematic  errors  in  the  measurement  of  £  discussed  above  are 
common  to  the  three  cases. 


The  coaperiton  between  data  and  the  model  predictiona  show  tba t  with  the 
addition  of  two  variables  to  describe  the  non-equil ibrina  kinetics  and  thermo¬ 
chemistry,  it  now  appears  that  finite-rate  radical  pool  kinetics  can  be 
predicted  with  reasonable  accuracy  in  the  contex*.  of  tnrbnlent  non-preaixed 
syngas  flaaes.  Results  plotted  here  show  better  agreement  with  (re-analyzed) 
data  on  aaj or  species  than  previously  thought.  Thus,  this  study  shows  that 
the  partial-equil  ibr iua  model  for  the  oxyhydrogen  radical  pool  including  CO  is 
more  useful  than  was  concluded  by  Correa  (1986a)  and  Pope  and  Correa  (1986). 
Although  there  are  unresolved  discrepancies  between  the  predictions  and  the 
data,  these  are  perhaps  due  to  the  difficulty  of  obtaining  and  interpreting 
Raman  data  to  a  greater  extent  than  previously  thought.  The  re-analyzed  data 
do  not  show  large  amounts  of  CO  present  where  is  absent,  which  would  have 
been  a  strong  indication  of  the  breakdown  of  partial  equilibrium.  Generally 
good  agreement  on  the  minor  species  (OH)  in  this  flame  obtained  earlier 
(Correa  et  al .  (1984)]  further  supports  the  partial-equilibrium  model  and 

emphasizes  the  importance  of  the  relatively  slow  rad leal -recombine tion  chemis- 


3.4  High  Reynold*  Nuaber  Flue:  Selection  of  Fuel 

One  of  the  pr inary  goal*  of  this  etndy  it  to  obtain  a  better  understan¬ 
ding  of  non-equil ibr ina  cheaiatry  effects  in  turbulent  jet  flaaet  at  high 
strain  rates.  Snch  flaaet  aay  have  regions  of  locally  extinguished  eddies 
eabedded  in  the  flaae,  which  aay  be  a  precursor  to  flaae  lift-off  and  blow- 
off.  The  probability  of  finding  large  regions  of  local  extinction  in  the 
flaae  it  a  strong  function  of  the  fuel  being  used.  The  choice  of  the  fuel  was 
priaarily  aotivated  by  the  desire  to  inorease  this  probability.  However, 
there  were  other  consideration  involved  in  the  choice  of  the  fuel,  which  are 
discussed  below.  The  fuel  chosen  for  detailed  study  consists  of  10%  ,  40% 

CO  and  30%  (by  voluae)  and  is  referred  to  as  LHC  (Low  Hydrogen  Content)  gas 
hereafter.  The  priae  considerations  in  the  selection  of  the  fuel  were: 

1.  The  probability  of  extinction  in  the  turbulent  flaae  should  be  suffi¬ 
ciently  high  to  warrant  detailed  aeasureaents.  To  estiaate  the  proba¬ 
bility  of  extinction,  a  siaple  experiaent  was  conducted  to  estiaate  the 

value  of  the  critical  strain  rate  paraaeter  (a  -2  V/a).  A  counter-flow 

c 

diffusion  flaae  burner  of  the  type  used  by  Tsuj 1  and  Yaaaoka  (1970)  and 

described  in  Drake  (1913)  was  used  to  estiaate  the  value  of  a  for 

c 

various  fuel  coapo sit  ions.  The  burner  consists  of  a  5  ca  diaaeter 
cylinder  (3  oa  long)  with  a  40*  sector  porous  plug  which  it  used  to  sta¬ 
bilise  a  flaae  2-3  aa  away  froa  the  surface  of  the  porous  burner.  The 
burner  it  installed  in  the  subsonic  low  turbulence  varlable-free-streaa 
velocity  tuanel  such  that  the  60  '  sector  is  syaaetric  aboat  the  stagna¬ 
tion  atreaa  line.  Aa  the  free-streaa  air  velocity  is  increased,  the 
flaae  it  pashed  closer  to  the  burner  surface  and  eventually  blows  off 


abruptly  at  a  critical  strain  rate,  for  which  the  characteristic  airing 
(iBe  evidently  exceeds  the  characteristic  reaction  tine  r c . 

Figure  9  shows  a  plot  of  ac  as  a  function  of  percent  Hydrogen  in 
the  fuel  which  consists  of  501  ^  and  the  rest  CO.  In  these 
aeasureaents  the  fuel  velocity  ratio  was  held  constant  to  avoid  exces¬ 
sive  heat  loss  froa  the  flaae  when  it  is  close  to  the  burner.  As  can  be 

seen  in  the  Figure,  the  value  of  a  increases  froa  400  s_1  corresponding 

c 

to  2%  to  1350  a*1  corresponding  to  15*  ^ .  The  value  of  oc  is  950 
s'1  for  the  LHC  gas  chosen  for  this  study.  This  value  is  auch  lower 
than  the  estiaate  (12,000  s-*)  for  Hydrogen  and  the  estiaate  (greater 
than  2000  s-1)  for  aed  BTC  gas  (40*  CO,  30*  H^,  30*  )  and  is  of  the 

order  of  oc  for  aethane  flaaes  (-400s  • 

The  probability  of  local  extinction  thus  decreases  very  rapidly  as 
the  hydrogen  content  is  increased.  For  pure  hydrogen  flaae  it  is  alaost 
iapossible  to  achieve  extinction.  Even  though  data  presented  by  Dibble 
and  Magre  (1987)  show  existence  of  eddies  that  are  lean  aixtures  of 
air-fuel  sad  have  very  low  teaperatur es,  these  eddies  aay  be  the  result 
of  large  holes  in  the  flaae  through  which  air  has  penetrated  into  the 

aain  jet  and  aixed  with  fuel  without  burning.  In  fact,  that  flaae  was 

visually  unsteady.  la  hydrocarbon  flaaes  there  is  evidence  of  local 
extiaetiou  at  relatively  low  strain  rates  (as  suggested  by  the  low  value 

of  «e)  but  hydrocarbon  flaaes  have  the  eoapl ications  of  Involving  com- 
plex  cheaistry  which  is  not  well  understood  and  of  being  extreaely 
aaaaeaable  to  Kaaan  aeasureaent  due  to  fluorescence  interferences  [Dib¬ 
ble  et  al .  (1987)].  These  aeasureaents  of  a  indicate  that  for  the  LHC 


gas  chosen  for  this  study  the  probability  of  extinction  is  fairly  signi- 
f  icant . 

The  fuel  should  have  relatively  simple  chemistry.  The  finite  rate  chem¬ 
istry  effects  of  flames  are  fairly  well  understood.  Both  H2 

and  H^/CO  turbulent  flames  at  low  Reynolds  numbers  have  been  studied  in 
detail  in  the  past  [Drake  and  Kollmann  (1985)].  The  study  of  H^/CO 
flames  is  also  important  because  most  hydrocarbon  fuels  undergo  pyro¬ 
lysis  and  form  «nd  CO  as  intermediate  constituents  before  burning. 

The  resulting  flame  should  be  amenable  to  Raman  diagnostics.  Even 
though  it  is  advantageous  to  reduce  the  content  of  the  fuel,  it  was 
found  that  for  less  than  10%  content,  the  amount  of  CO  had  to  be 
increased  significantly  to  obtain  a  stable  flame.  Excessive  amounts  of 
CO  present  in  the  flame  seriously  interfere  with  Raman  diagnostics  due 
to  chem i-1  urn ine scent  interference  from  GO+O^CO^  reaction.  This  chemi¬ 
luminescence  is  fairly  broad-banded  and  extends  all  the  way  to  the  far 
red.  Even  though  flame  luminescence  can  be  subtracted  on  a  shot-to-shot 
basis  during  Raman  data  collection,  the  change  in  the  chemi-luminescence 
during  the  time  interval  can  result  in  spurious  Raman  signals.  Thus,  it 
was  necessary  to  reduce  both  the  CO  end  ij  content  of  the  flame  as  much 
as  possible  by  diluting  it  with  Nitrogen.  For  Nitrogen  dilution  more 
than  50%  a  stable  flame  at  high  Reynolds  number  could  not  be  obtained 
even  with  the  aid  of  a  pilot.  The  LHC  flame  at  15,000  Reynolds  number 
was  fairly  stable  with  the  premixed  pilot.  Thus,  the  Raman  feasibility, 
the  stability  of  the  flame  and  the  probability  of  localised  extinction 
required  a  compromise  in  the  choice  of  the  fuel  composition. 
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A  final  consideration  regarding  the  choice  of  the  fuel  va s'  the  ability 
to  nse  Rayleigh  scattering  to  obtain  2D  images  of  density  in  the  flame 
by  matching  the  Rayleigh  cross-section  of  air,  fnel  and  products.  With 
a  fuel  composed  of  10b  t  40%  CO,  25b  Nj,  and  25b  Ar  it  is  estimated 
that  the  change  in  Rayleigh  ross-section  between  the  fnel,  air  and  pro¬ 
ducts  is  less  than  7b.  This  is  an  important  property  of  the  modified 
LHC  gas  for  the  measurement  of  strain  rate  in  the  future. 
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3.5  Pilot-Stabilized  Burner 


The  batic  subsonic  1 ov- turbul ence  combustion  tunnel  has  been  described  in 
detail  elsewhere  (Lapp  et  al .  (1983)).  The  combustor  has  been  used  exten¬ 
sively  for  experimental  studies  of  laminar,  transitional  and  turbulent  Hydro¬ 
gen  diffusion  flames  in  the  past.  The  combustor  was  redesigned  and  rebuilt  to 
stabilize  the  high  Reynolds  number  LHC  flame  by  providing  a  coaxial  premixed 
pilot.  The  premixed  pilot  burner  is  designed  in  a  manner  similar  to  that  sug¬ 
gested  by  Sterner  and  Bilger  (1985).  The  turbulent  jet  diffusion  flame  is 
stabilized  at  the  tip  of  a  nozzle  of  diameter  3.18  mm  by  the  co-annular 
premixed  pilot  flame.  The  flow  rates  are  measured  using  critical  flow  orifi¬ 
ces.  Two  flames  at  cold  flow  exit  velocities  of  25  m/s  and  80  m/s  correspon¬ 
ding  to  Reynolds  numbers  of  5000  and  15,000,  respectively,  were  studied. 

The  high  Reynolds  number  turbulent  flame  of  LHC  fuel  (ReK15,000)  provides 
a  fluctuating  strain  field  which  can  cause  local  extinction  when  some  critical 
strain  rate  is  exceeded  locally.  Extensive  regions  of  such  localized  extinc¬ 
tion  can  result  in  eventual  liftoff  or  blowoff  of  the  flame.  The  pilot- 
stabilized  burner  is  designed  to  provide  a  stable  flame  in  the  region  of  large 
strain  near  the  exit  yet  devoid  of  gross  unsteady  features.  The  pilot  flame 
results  in  reduced  initial  strain  rates  and  provides  radicals  resulting  in  the 
stabilization  of  the  main  flame.  Further  downstream  the  turbulence  is  fully 
developed  and  the  strain  rate  is  still  high.  Figure  10  shows  a  schematic  of 
the  burner  exit.  The  pilot  flame  consists  of  a  stoichiometric  mixture  of  the 
fnel  ( 10%  Hj,  40%  CO,  501  Nj)  and  air  which  were  mixed  well  upstream 
of  the  combustor  exit.  This  fuel-air  mixture  ratio  (ts«0.43)  results  in  the 
combustion  products  of  the  premixed  flame  having  the  same  C/H/0  ratios  as  the 


main  fuel  when  bnrned  stoichiometr ical ly  and  simplifies  the  modeling  of  the 
flowfield  by  retaining  the  two-stream  feature.  The  premixed  annular  burner 
has  an  open  area  five  times  that  of  the  main  jet.  A  flame  arrestor  consisting 
of  18  holes  of  1  mm  diameter  each  is  installed  4  mm  upstream  of  the  exit  to 
anchor  the  premixed  flame.  The  exit  velocity  of  the  premixed  products  is 
estimated  to  be  20  mj  s  based  on  an  estimated  post-flame  temperature  of  2000*K. 
The  co-flowing  air  velocity  in  the  tunnel  was  set  at  5  m/s.  Figure  11  shows  a 
schematic  of  the  flame  identifying  the  various  regions  in  the  flame.  The 
pilot  flame,  the  main  jet  and  the  region  of  maximum  strain  rate  where  the 
flame  is  narrow  are  also  identified  in  the  figure. 

The  main  jet  flame  at  Re  *=  15,000  is  approximately  50  diameters  long 
(based  on  main  jet  diameter).  The  premixed  pilot  flame  extends  to  about  5-6 
diameters  when  lit  without  the  main  jet.  It  is  estimaMl  that  the  enthalpy 
content  contributed  by  the  pilot  is  less  than  5%  of  that  bj  the  main  jet  due 
to  its  relatively  low  volumetric  flow  rate. 


3.6  Improvements  in  Raman  System 

Dae  to  the  high  CO  content  of  the  flame  and  the  high  Reynolds  number,  it 
*as  found  that  there  was  severe  interference  from  chemi-luminescence  of  the 
flame  (blue)  in  all  channels  of  the  Raman  system  including  N^(AS),  {^(S),  O2. 

1^0  and  H2 .  The  background  flame-luminescence  was  as  large  as  tbe 
associated  Raman  signals  in  most  cases.  For  the  Raman  diagnostic  system  to 
acquire  meaningful  data,  not  only  should  the  actual  background  levels  be  low, 
but  the  change  in  flame  background  over  the  integration  period  of  the  Raman 
signal  should  be  low  compared  to  the  signal  itself.  The  following  major 
improvements  were  made  in  the  Raman  system  to  minimize  the  problems  due  to 
interference  from  flame  luminescence. 

1.  Introduced  a  polarization  filter  in  the  collection  optics.  Since  the 
Raman  scattered  light  is  polarized  (98%)  at  the  same  polarization  as  the 
incident  laser  light,  whereas  the  flame  chemiluminescence  is  unpolar¬ 
ized,  a  factor-of- two- improvement  in  signal-to-noise  ratio  was  achieved 
by  the  introduction  of  the  polarization  filter. 

2.  The  laser  was  focussed  tightly  by  using  a  smaller  focal  length  lens. 
This  enabled  a  reduction  in  the  entrance  slit  size  and  a  corresponding 
reduction  in  the  flame  chemi-lumi.uescence  without  reducing  the  Raman 
signal . 

3.  By  reducing  the  elapsed  time  between  the  collection  of  data  end  the 
firing  of  laser  beam  from  13  psec  to  8  psec,  the  effect  of  change  in 
flame  background  daring  Raman  integration  time  on  the  Raman  signal  was 


minimized. 
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4.  A  new  photomrl t ipl ier  tube  to  monitor  the  Rayleigh  scattering  signal  was 
introduced  in  the  systeai.  Since  the  N2  channel  is  affected  most  by  the 
luminosity  of  the  flame,  the  temperature  measurements  using  the  Stokes- 
Anti-Stokes  method  are  susceptible.  Bence.  Rayleigh  scattering  was  used 
to  obtain  an  independent  measurement  of  temperature.  The  Rayleigh  scat¬ 
tering  intensity  Ig  i*  given  by  the  following  expression  (Dibble  and 
Hollenbach  (1981)). 

I,  =  CINQ1  (Jg)  (2) 

where,  C  is  a  calibration  constant  of  the  collection  optics,  1  is  the  incident 
laser  intensity,  N  the  molecular  number  density,  0  the  solid  angle  of  the  col¬ 
lection  optics,  a  is  the  Rayleigh  scattering  setting  cross-section  of  the  gas 
molecules  and  1  the  length  of  the  laser  beam  segment  imaged  onto  the  detector. 
In  a  uniform  pressure  field,  the  above  expression  is  simplified  using  the 
ideal  gas  law  to  the  following  form: 

J,  “  CIQ1  If  <ft>  “  T  ^dQ^ef f.  (3) 

where  I  is  a  constant,  P  is  the  pressure,  T  is  the  temperature  of  the  gas  and 
R  is  the  universal  gas  constant.  The  above  expression  shows  that  variations 


in  Rayleigh  scattering  intensity  result  both  from  temperature  variations  and 
changes  in  scattering  cross-sect  ion.  However, 


(d S.) 

dQ'eff 


and  hence  1^  ■  which  implies  that  the  Rayleigh  scattered  intensity  can 
be  related  to  temperature  directly  in  jet  diffusion  flames  only  if  either 
additional  information  about  the  change  in  the  Rayleigh  cross-section  of  the 
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gas  mixture  is  available  or  if  the  average  Rayleigh  cross-section  is  constant. 
In  this  flame  the  Raman  scattering  signals  from  the  major  species  were  used  to 
obtain  information  about  changes  in  the  Rayleigh  cross-section  of  the  mixture 
and  thus  measure  temperature  in  an  iterative  manner. 


3.7  Experimental  Results 


Raman  measurements  of  instantaneous  temperature,  mixture  fraction  and 
mole  fractious  of  major  species  were  obtained  in  two  flames  at  Re  =  3,000  and 
for  the  Re  -  13,000.  No  pilot  flame  was  used  for  Re  *  3,000  flame.  The  tem¬ 
perature  was  measured  using  three  methods  (1)  St okes-Anti- Stokes  method  based 
on  nitrogen,  (2)  sum  of  mole  fractions  of  major  species  and  (3)  Rayleigh  scat¬ 
tering  as  discussed  earlier.  Figure  12  shows  a  comparison  between  the  mean 
centerline  temperature  in  the  Re  *  15,000  flame  obtained  using  the  above  three 
methods.  The  variation  in  the  mean  temperature  obtained  using  the  three 
methods  is  at  most  100*K.  The  Rayleigh  scattering  channel  corresponds  to  the 
best  signal-to-noise  ratio  and  is  the  most  reliable  measurement.  The  mixture 
fraction  was  obtained  based  on  the  carbon  element  using  the  following  defini¬ 
tion. 


XCO  +  XCO„ 


5  “  <- 


(*iwi) 


major  species 


Figures  13,  14  and  15  show  the  scattergrams  of  temperature  vs.  mixture 
fraotion  at  various  locations  in  the  two  flames.  Each  data  point  in  the 
figures  represents  an  instantaneous  simultaneous  realization  of  temperature 
and  mixture  fraction  (based  on  carbon  element).  Figure  13  is  the  scattergram 
at  x/d-10,  r/a-1.6  in  the  Re-3000  LHC  flame  ("a”  refers  to  the  radius  of  the 


■sin  jet;  a  =  1.59  bib).  The  solid  line  in  the  figure  represent s  an  adiabatic 
equilibria®  (AE)  calculation.  As  can  be  seen  in  the  Figure,  aost  of  the  data 
lies  on  the  slightly  fuel  rich  side  of  the  curve  and  the  scatter  is  ainiaa]  . 
The  teaperature  variation  froa  AE  is  less  than  500*K  at  this  location.  The 
departures  froa  equilibria®  represent  the  effect  of  turbulence  and  aeasureaent 
error.  However,  the  departures  froa  equilibria®  are  saal 1  and  do  not  present 
evidence  of  extinction. 

Figures  14  and  15  present  data  obtained  in  the  high  Re*=15,000  flaae  at 
x/d=10,  r/ a  =  1.89  and  x/da60,  r/a  =  3.77  locations  respectively.  These  loca¬ 
tions  are  representative  of  the  regions  of  high  strain  rate  and  of 
coapletely-aixed  equilibria®  flow,  respectively.  Figure  14  shows  an  enoraous 
aaount  of  scatter  in  the  data  (2000  points)  and  large  departures  froa  equili¬ 
bria®  values.  The  airfare  fractions  lie  in  the  range  0.3  to  0.6  and  teapera¬ 
ture  departures  froa  equilibria®  in  the  range  400-1200*K.  The  data  points 
having  aixture  fraction  values  near  stoichioaetr ic  and  significantly  low  tea- 
peratuxes  aay  represent  locally  extinguished  flaaelets.  Local  extinction 
refers  to  a  stoichioaetr  ic  aixture  of  fuel  and  air  which  is  unburnt  or  has 
been  extinguished;  however,  the  data  indicates  liaited  evidence  of  local 
extinction  i.e.,  in  approximately  only  20%  of  the  realizations.  Figure  15 
shows  the  scattergraa  at  x/d**40  at  which  location  the  fuel  and  air  are  signi¬ 
ficantly  aixed  and  the  strain  rates  are  significantly  lower.  As  can  be  seen, 
aost  of  the  data  lies  along  the  fuel  rich  side  of  the  AE  curve  and  the  depar¬ 
tures  froa  equilibria®  are  minimal  even  though  there  is  significant  scatter  in 
the  data.  There  is  virtually  no  evidence  at  all  of  localised  extinction  at 


this  axial  location. 
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3.8  Lu i  nx  r  Flaaelet  Model:  Predictions  and  Comparison 

Under  the  assumption  that  the  laainar  flaae  thickness  is  saal 1  compared 
with  the  tnrbnlent  (Kolmogorov.  Batchelor)  scales  in  the  flow,  the  "laainar 
flaaelet"  analysis  for  non-preaixed  turbulent  flames  has  been  advanced,  e.g.  , 
by  Lie*  et  al .  (1981).  In  this  aodel.  the  Instantaneous  composition  of  the 
turbulent  flow  is  assuaed  to  correspond  to  that  in  an  undisturbed  laainar  dif¬ 
fusion  flaae  at  the  prevailing  value  of  mixture  fraction.  As  the  authors 
pointed  out,  this  aodel  is  superior  to  chemical  equilibrium  models,  particu¬ 
larly  in  fuel-rich  hydrocarbon  flaaes  where  CO  levels  are  far  froa  equili¬ 
brium.  The  flaaelet  aodel  also  offers  an  alternative  to  explicit  aulti- 
variable  representations  of  the  chemistry,  such  as  performed  by  Janicka  and 
Kollaan  (1979)  and  Bilger  (1980)  for  turbulent  diffusion  flaaes  and  Correa 
et  al.  (1984)  for  00/ flames.  The  undisturbed  flaaelet  aodel  has  also  been 
used  in  other  analyses,  eg.,  Faeth  and  Saaul son  (1983)  for  sooting  flaaes  and 
Bazdan  and  Stevens  (1983)  for  CO  flaaes.  Turbulent  flaaes  are  hypothesized  to 
consist  of  an  ensemble  of  these  laainar  flaaelets. 

Numerical  studies  of  the  structure  of  laainar  diffusion  flaaes  with  fill 
kinetics  indicate  a  strong  dependence  on  the  strain-rate  (Miller  et  al. 
1984).  A  configuration  often  studied  is  the  counterflow  porous-cylinder  bur¬ 
ner  as  used  by  Tauj i  and  Taaaoka  (1970).  It  is  now  established  that  as  the 
strain  rate  2V/ a  -  where  "V"  is  the  oxidizer  velocity  and  "a"  the  cylinder 
radius  -  is  increased,  the  aaxiaua  teaperature  in  the  system  decreases  until 
an  abrupt  collapse  to  inert  flow  occurs,  at  soae  critical  level  of  strain.  The 
critical  strain  varies  with  fuel  type  ranging  froa  ~  12,000  s~*  for  [David 
et  al .  (1987)]  to  *>400  s~*  for  aethane  (Miller  et  al.  1984),  and  is 


approximately  850  s  1  for  the  LUC  ga  s  chosen  for  this  study. 


The  usportsnce  of  the  scslsr  dissipstion  rste  to  flame  structure  csn  be 
shown  ss  follows.  When  the  species  conservation  equation 


dY.  .  dY, 

pu  ~  “  «  (D.  — l)  +  ir . 

k  ai^  1  dx^  i 


(5) 


is  transformed  to  a  mixture- fract ion  space,  i.e..  a  coordinate  perpendicular 
to  surfaces  of  constant  (elemental)  composition,  it  becomes 


D. 
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(6) 


which  shows  that  when  the  scalar  dissipation  x  -  Di  (:^£-)2  gets  too  large,  the 

‘2\ 

term  denoting  reaction  -  is  reduced  and  the  flame  can  be  extinguished. 

«u2 


The  sensitivity  of  the  laminar  flamelet  to  atrain  implies  that  the 
instantaneous  composition  in  a  turbulent  flame  will  depend  on  the  local  strain 
as  well  as  mixture  fraction.  This  has  been  accounted  for  by  Liew  et  al . 
(1984),  who  showed  that  the  joint  pdf  would  govern  the  instantaneous 

) 

structure.  In  physical  terms,  the  composition  would  correspond  to  that  at  the 
same  mixture  fraction  in  a  laminar  flame  subjected  to  the  same  strain.  In 

principle,  a  "library"  of  such  strained  flames  could  be  generated  computa- 

I 


tionally  or  experimentally;  however,  Liew  et  al .  (1984)  adopted  the  simpler 
approach  of  sampling  the  undisturbed  flame  (from  experimental  data)  or  inert 
flow,  depending  on  X  at  the  location  of  maximum  temperature  (X^  )  being  less 

or  more  than  an  assumed  quenching  value  of  scalar  dissipation,  1^,  In  highly 
turbulent  flames  where  the  local  pdf  of  strain  results  in  1>  1^  much  of  the 
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tiae,  the  kid  f  1  crw  would  then  be  strongly  biased  towards  the  inert  flow. 
Local  extinction  would  therefore  be  very  significant. 


The  aodel  used  in  this  study  follows  Liew  et.  al .  (1984)  but  assnaet  tbai 
equilibria*  (frosen)  flow  exists  instantaneously  depending  on  whether  the 
local  scalar  dissipation  rate  is  less  (nore)  than  soae  critical  value  X^.  lbe 
value  of  1^  u  critical  to  the  aodel  predictions  and  here  the  results  are  cob 
pared  for  a  range  of  values  of  X^. 

The  argaaent  by  which  the  local  distribution  of  scalar  dissipation  rate 
is  obtained  follows  Liew  et  al .  <1984).  Kolaogorov’s  hypothesis  of  a  log- 

noraal  distribution  is  aade  with  the  aean  Z  given  by 

z  -  c  r*'2  r/z  <7) 

h 

which  is  the  aodel  led  tera  in  the  transport  equation  for  the  variance  of  nix- 
tare  fraction.  The  locel  pdf  for  1  is 


P(X)  ■*  i  - i-jr-  exp[~(  lnI-|i)2/2o2] 

1  9  <2n>* 


where  o2  is  the  variance  of  InX.  The  pdf  of  X  conditioned  on  the  presence  of 
xones  of  intense  reaction  is  assuaed  to  be  identical  and  aay  be 

integrated  to  give  the  probability  of  being  below  a  critical  level  X^( 
the  probability  of  coabnation  P 


pc  -  i-i  *rfc  0" ®] 


(9) 


As  the  gas  flows  downstreaa,  the  strain  rates  decrease  and  the  probability 
of  the  gas  being  below  the  burning  liait  increases.  Thns.  wore  of  the  gas  is 


*q  ui  1  1 br  a  t  *d . 


The  pdf  over  theraocbea  ical  states  it  needed  to  obtain  various  qnentitici 
tocb  et  Mean  density.  Since  the  floe  it  either  inert  (o)  or  in  a  state  of 
equilibria*  tel.  the  pdf  ia  P(£)  -  ll-Pc)  Pu ( 4 )  ♦  P£  •  ( 4 )  where  P£  it  the 
fraction  of  gat  which  it  burning.  The  pdf  for  the  occorrence  of  valoet  of  the 
■  lxtnre  fraction  it  needed.  Soch  pdf  *  a  can  be  obtained  directly  fro*  evolu¬ 
tion  equations  bated  on  the  Nav le r— St oke a  and  apeciea  equations  [Pope  (1981)1 
or  more  tiaply  by  assuaed-thape  pdf  Bethods  at  in  Kent  and  Bilger  (1976)  for 
P(4). 

The  pdf  it  widely  edopted  at  being  a  be ta- funct ion 

P<  4)  -  p(()  ( 10> 

which  can  be  pa  rue  ter  it  ed  in  terui  of  ita  lowest  aoaents,  i.e.,  the  scan  and 
variance,  and  ao  varies  throughout  the  flowfield. 

The  coapnta tional  procedure  aay  be  snaaerised  as  follows:  The  shear- 

layer  equations 

ii  (pu#i>  4  r  h  (rp^i>  ■  i  h  uri  af)  *  \  (11) 

are  solved  in  x-etreaa  function  coordinates  nsing  the  nuaerical  technique  of 
Patankar  and  Spalding  (1972).  Here  9^  is,  respectively,  the  aiial  velocity, 
turbulence  kinetic  energy,  dissipation  rate,  and  Bean  and  variance  of  aixture 
fraction.  The  source  t«rt  and  the  transport  coefficient  are  chosen 
appropriately  for  particular  equation.  Inlet  and  boundary  conditions  on  the 
variables  correspond  to  fully-developed  turbulent  fuel  flow  with  an  assuaed 
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intensity  end  length-scale  of  turbulence.  The  flow  exiting  the  annulus  is 
essaaed  to  be  coapletely  burned. 

At  each  point  in  the  floe  the  local  properties  are  obtained  by  weighting 
between  the  inert  and  burning  flow  properties  with  the  probability  of  coabus- 


tion,  P 

c  ’ 

obtained  froa  Eqs. 

(d).  The 

properties  of 

the  burning 

are 

obtained 

by 

convoluting  the 

local  pdf 

PU,?) 

froa 

Eqs.  (10) 

with 

the 

corresponding 

theraochea  ical 

properties. 

such  as 

KO 

or  p((,n). 

froa 

the 

equilibriua  coabustion  aodel  . 

Figure  16  shows  the  centerline  axial  variation  of  aean  teaperatnre.  The 
aodel  predictions  for  «  40,  100  and  1000  s_1  ere  also  shown.  The  exper- 

iaental  data  follows  a  siailar  trend  though  the  aodel  predicts  a  slower  growth 

rate  for  all  values  of  X^,  The  deviations  of  predictions  froa  data  are 
greatest  at  x/d  "  20.  Further  downstreaa  (e.g.,  at  x/d  *  30)  the  strain  rates 
are  expected  to  be  auch  lower  and  so  the  predicted  profiles  for  all  these 

values  of  X^  collapse  and  arc  auch  closer  to  the  experiaental  results.  This 
however,  is  not  conclusive  support  of  the  aodel  because  if  a  flaaelet  is  quen¬ 
ched  in  the  large  strain  region  (lower  x/d  values)  then  it  is  preaixed  and 
cannot  be  described  adequately  by  the  aodel.  This  is  because  the  aodel  is 
applicable  to  strained  flaaes  of  the  diffusion  type  only  whereas  extinguished 
flaaelets  lead  to  a  preaixed  situation. 

Radial  profiles  of  the  aodel led  scalar  dissipation  rate  for  X^  -  40,  100 
and  1000  at  the  axial  location  x/d  ■  20  (Figure  17)  ahow  large  differences. 
This  is  because  the  jet  is  closer  to  equilibriua  at  the  higher  1^  which 

changes  the  penetration  and  aixlng.  These  differences  are  reflected  in  signi- 
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ficant  variations  of  the  probability  of  conbustion,  i.e.,  the  expectation  1<X 
(Figure  18).  Although  it  is  not  at  present  possible  to  aeasure  the  dissipa¬ 
tion  rate  -  which  str ic t ly- speaki ng  would  require  Measuring  the  aixture  frac¬ 
tion  in  three  diaensions  with  sufficient  spatial  and  teaporal  resolution  to 
resolve  the  aaallest  eddies  -  the  scattergraas  of  teaperature  vs.  Mixture 
fraction  in  the  zone  of  greatest  strain  indicate  approxiaately  a  20%  probabil¬ 
ity  of  local  extinction,  as  discussed  above  (Figure  15).  This  is  of  the  same 
Magnitude  as  the  probabilities  in  Figure  18  for  1^  .  100.  The  region  of  aax- 
iaua  probability  of  extinction  and  the  trand  of  its  radial  variation  are 
predicted  well  by  the  aodel . 

Figures  19-21  show  the  density-weighted  Bean  aixture  fraction,  conven¬ 
tional  Beans  of  density  and  teaperature  at  x/d“20.  The  jet  spreads  less 
rapidly  than  predicted  by  the  aodel  as  seen  by  the  narrower  profile  of 
aixt  ur  e-  f  ract  ion  (Figure  19}  and  the  centerline  profile  of  aixture  fraction 
(not  shown).  Differences  in  these  profiles  for  different  X^  are  partially 
responsible  for  the  differences  in  the  aean  density  and  teaperature  (Figures 
20  and  21),  but  the  latter  differences  are  also  due  to  the  differing  dissipa¬ 
tion  rates  and  probabilities  of  coabustion.  Thus  the  choice  of  X^  has  a  large 
influence  on  the  results.  Froa  Figure  19  it  appears  that  the  best  coaparison 

is  obtained  for  I^»1000,  but  that  iapl  ies  an  equilibriua  flaae  and  is  in  con¬ 
tradiction  with  the  observations  aade  in  Figures  16  and  18.  The  agreeaent  for 

1^*1000  should  be  viewed  as  fortuitous  and  probably  due  to  errors  in  the 
spread-rate  due  perhaps  to  ais-est iaates  of  inlet  conditions.  Veil  downstreaa 
of  the  flaae  xone,  the  strain-rate  is  reduced  and  predictions  becoae  rel  a- 
insensitive  to  X  Figures  22-24  again  coapare  predicted  radial  pro- 


filet  of  density-weighted  tern  mixture  fraction  and  conventionally-averaged 
mean  density  and  temperature  with  the  Raman  data.  The  probabilities  of  com¬ 
bustion  tend  to  unity  in  all  the  cases  -  because  the  strain-rates  are  so  low 
at  this  location.  The  predicted  profiles  are  quite  close  to  each  other  and  to 
the  data.  The  relatively  good  agreement  with  the  data  may  be  fortuitous 
because  as  discussed  earlier,  the  model  does  not  represent  premixed  laminar 
flamelets  (which  were  earlier  quenched  in  regions  of  high  strain  rates)  ade¬ 
quately. 

Discrepancies  in  the  region  of  greatest  strain  near  the  nozzle  exit  may 
be  due  to  one  or  more  of  several  possible  causes.  Exit  conditions  in  the  noz¬ 
zle  and  pilot  are  difficult  to  characterize  and  may  exert  a  non-trivial 
influence  on  the  near-exit  region  despite  the  large  shear.  These  conditions 
include  turbulence  quantities  and  the  heat  loss  in  the  pilot.  From  model  cal¬ 
culations  it  was  found  that  the  results  are  sensitive  to  the  temperature  pro¬ 
file  at  the  exit  of  the  premixed  pilot  burner.  At  present  experiments  are 
being  conducted  to  measure  the  temperature  profiles  at  the  exit  of  the  pilot 
burner  directly  to  provide  the  initial  conditions. 

These  findings  may  be  summarized  as: 

1.  There  is  no  firm  evidence  of  local  extinction  (to  frozen  conditions)  in 
the  LHC  gas  Re  ■  15,000  flame  even  in  the  regions  of  largest  strain. 
The  model  predicts  the  maximum  probability  of  extinction  to  be  30%  in 
this  flame  at  x/d-20.  Experimentally,  it  is  observed  that  approximately 
20%  of  data  points  have  mixture  fraction  values  close  to  stoichiometric 
but  have  significantly  lower  temperatures  than  adiabatic  equilibrium. 
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These  low  temperature  eddies  could  result  from  any  of  the  following 
reasons:  (a)  they  could  represent  eddies  which  were  extinguished  but  are 
above  300*K  due  to  nixing  with  surrounding  hot  products,  (b)  they  could 
represent  strained  eddies  with  such  lower  temperatures  (than  AE)  due  to 
large  non-equilibrium  chemistry  effects,  (c)  they  could  represent  mixed 
eddies  which  are  partially  burned  (i.e.  they  were  extinguished  due  to 
large  strains  after  initial  ignition),  (d)  they  could  represent 
artifacts  of  Raman  probe  averaging  (Raman  probe  volume  is  200p  x  200p  x 
650p) .  Further  work  is  needed  to  identify  the  exact  nature  of  these 
data. 

2.  The  theoretical  predictions  of  the  model  show  that  the  results  are  very 
sensitive  to  the  assumed  value  of  X^.  the  quenching  parameter,  especi¬ 
ally  in  the  region  of  large  strain  and  are  also  sensitive  to  the  assumed 
velocity  and  temperature  profiles  at  the  exit  of  the  pilot  burner.  In 
the  next  year  of  the  program  the  exit  velocity  and  temperature  profiles 
near  the  tip  of  the  combustor  will  be  measured  using  laser  velocimetry 
and  Raman  scattering. 

3.  Even  though  the  model  provides  significant  insight  into  the  quenching 

and  localised  extinction  phenomena  in  turbulent  diffusion  flames,  there 
is  significant  disagreement  between  the  data  and  calculations  especially 
at  x/d-20.  The  relatively  better  agreement  at  x/d«50  may  be  fortuitous. 
The  model  suffers  from  serious  drawbacks  of  (a)  relying  on  a  relatively 
arbitrary  quenching  parameter  X^,  (b)  an  inability  to  adequately 

describe  flamelets  after  they  are  quenched  (premixed  strained 
flamelets),  and  (c)  considering  only  two  types  of  flamelets)  equilibrium 
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3.9  Alternate  Approach  to  Modeling  Turbulent  Extinction 

An  important  shortcoming  of  the  laminar  flamelet  analytes  is  that  in 
practice  laminar  flames  are  not  necessarily  thin  compared  with  turbulent 

scales.  In  and  CO/H2  flames,  there  are  two  distinct  sets  of  scales  associ¬ 
ated  with  the  chemistry,  that  corresponding  to  two-body  chain-branching  or 
propagating  reactions  and  that  corresponding  to  three-body  recombination  reac¬ 
tions.  The  large  disparity  between  these  scales  has  lead  to  partial- 
equilibrium  models  which  assume  the  two-body  reactions  to  be  in  equilibrium 
[Janicka  and  Kollman  (1979),  Bilger  (1980),  Correa  et.  al .  (1984),  Correa 

(1986a),  Pope  and  Correa  (1986)].  This  disparity  also  implies  that  in  intense 
turbulence  the  two-body  reaction  zones  would  remain  comparatively  thin  while 
the  three-body  zones  would  overlap,  leading  to  distributed  reaction  zones 
which  are  incompatible  with  the  laminar  flamelet  hypotheses. 

Thus,  if  independent  control  over  the  turbulent  mixing  was  exercised,  as 
the  time-scale  is  progressively  reduced  from  large  values  the  OB  (to  pick  a 
radical)  would  first  increase  as  the  recombina tion  reactions  are  no  longer 
fast  enough  to  realize  equilibrium,  and  then  abruptly  drop  as  the  two-body 
reactions  are  quenched  and  the  gas  is  inerted.  Temperature  shows  a  monotonic 
decrease  throughout  this  process  with  a  more  rapid  decrease  as  extinction  is 
approached.  This  behavior  is  consistent  with  the  notion  that  the  abrupt  tran¬ 
sitions  in  a  stretched  flamelet  are  characteristic  of  high-activation  energy 
(shuffle)  reactions  while  well-distributed  reaction  zones  are  characteristic 
of  low  or  zero  activation-energy  (recombination)  reactions. 

The  large  difference  in  the  rates  of  the  two-body  and  the  three-body 


reactions  is  apparent  from  experiments  at  relatively  low  Reynolds  numbers. 
Correa  et  al.  (1984)  nsed  the  partial-equilibrium  model  for  00/ in  an  Re  = 
8500  (cold  jet)  flame.  The  data  indicated  negligible  or  no  local  extinction  - 
i.e.  ,  none  of  the  Raman  data  shoved  lov  temperature  in  flammable  gas  -  but 
large  amounts  of  superequil  ibr  ium  OH,  up  to  factors  of  5.  The  model  agreed 
veil  vith  these  data.  One  may  conclude  that  the  time-scales  of  turbulence  in 
such  flames  are  short  enough  to  interfere  vith  three-body  recombination  chem¬ 
istry,  i.e.  a  Damkohler  number  based  on  the  recombination  reactions  is  of 
0(1).  The  time-scales  are  still  too  long  to  disturb  the  tv o— body  reactions 
^2-body^^^  lnd  *°  there  is  little  or  no  local  extinction. 

Here  it  is  proposed  to  retain  the  partial  equilibrium  model  but  vith  the 
quenched  flamelet  approach  to  the  tvo-body  reaction  zone.  The  tvo-body  reac¬ 
tion  zone  is  assumed  to  be  either  in  equilibrium  or  quenched  (i.e.,  leading  to 
inert  flov)  depending  on  the  local  strain  X  being  less  or  more  than  the  quench 
level,  X^. 

The  kinetic  mechanism  adopted  consists  of  five  fast  shnffle  reactions 
(S1)-(S5),  vhich  are  considered  to  be  in  partial  equilibrium. 


H  +  °2  *  OH  +  0  (SI) 

0  +  ^  =  OH  +  H  ( S2) 

^  +  OH  -  HjO  +  H  (S3) 

20H  *1^0  +0  (S4) 

H  +  H02  .  OH  +  OH  ( S5) 


and  the  four  relatively  slov  three-body  recombination  reactions  (Rl)  -  (R4) 
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H  +  OH  +  M  — >  H^o  +  M  (Rl) 

H+O+M— >OH+M  (R2) 

H+H+M-^^  +  m  (R3 ) 

H  +  02  +  u  ->  H02  +  M  (R4) 


First  it  mast  be  recognized  that  the  scalar  dissipation  rate  continuously 
decreases  as  the  fluid  flows  downstream.  The  region  of  intense  stretch,  which 
is  visualized  as  being  co-located  with  the  2-body  shuffle  reaction  zone, 
therefore  experiences  a  transition  through  the  critical  strain  rate,  and  the 
2-body  reactions  are  assumed  to  iaaediately  reach  partial  equil ibr ium.  This 
is  analogous  to  the  approach  used  by  Liew  et  al.  (1984)  who  assuaed  a  burnt 
flaaelet  state  for  the  flaaelets  which  were  below  a  critical  strain  rate.  In 
the  present  case  only  the  2-body  reactions  are  assuaed  to  reach  equilibrium 
the  recoabination  reactions  coaaence  and  drive  the  system  towards  full  chemi¬ 
cal  equil ibr iua. 

The  concept  is  easily  visualized  as  follows  in  one  diaension.  If  the 
pdf's  for  scalar  disaipation  rate  X  at  an  upstream  and  downstream  location  are 
denoted  by  P&(X)  and  Pd(X),  respectively,  and  the  probabilities  of  2-body 
reaction  equilibration  are  *u(X<Xcr)  and  Pd  (X<Xcr),  then  the  probability  that 
the  recoabination  reactions  "just''  started  is  Pd  (I<Xcr)  -  PB  (I<Xcr).  The 
concept  is  currently  being  extended  to  aul ti-diaensional  flow  and  will  be  used 
to  replace  the  analysis  of  Section  3.8. 
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3.10  Future  Work 

For  the  next  year  of  the  program  the  following  tacks  have  been  identi¬ 
fied: 

1.  Characterize  the  temperature  and  velocity  profiles  at  the  exit  of  the 
pilot  combustor  using  Raman  scattering  and  laser  velocimetry. 

2.  Attempt  to  increase  the  probability  of  local  extinction  in  the  flame  by 
further  diluting  the  fuel  with  nitrogen  and  increasing  the  Reynolds  num¬ 
ber  of  the  flame  as  much  as  possible. 

3.  Make  Raman  measurements  in  the  Hj/CO  flame  to  help  identify  the  origin 
of  the  low  temperature  fuels-air  eddies  observed  in  the  flame. 

4.  Obtain  planar  OH/Rayleigh  images  in  the  flame  to  determine  flame 
connectivity/ strain  field  at  high  Reynolds  numbers. 

3.  Complete  the  modeling  of  the  high  Reynolds  number  plane  using  the  joint 
partial  equil ibr ium/ f lamel et  type  model  discussed  above  and  compare  with 


experimental  results 
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Lnt  of  Tablet 


1.  High  temperature  correction  factors  for  vibrational  Raman  scattering 
naed  to  correct  for  fraction  of  hot  bends  of  variona  species  falling 
outside  their  respective  exit  slits. 

2a.  Temperature  and  mole  fraction  of  major  species  at  x/d=25  obtained 
without  correcting  for  high  temperature  effects  (correction  factors 
set  to  unity). 

2b.  Temperature  and  mole  fraction  of  major  species  at  x/d“25  corrected 
for  high  temperature  vibrational  bends  nsing  the  correction  factors 
1 ist ed  in  Tabl e  1 . 

3.  Mean  mixture  fraction  at  y=0  obtained  based  on  the  Carbon,  Hydrogen  and 


Nitrogen  elements. 
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Fig.  18 
Fig.  19 
Fig.  20 


Axial  variation  of  neas  temperature  at  y  =  0  in  the' Be  *  15,000 
LHC  flame.  Model  predictions  for  I  0f  40,100  and  1000*1 
also  shown  in  the  figure.  4 

Badial  profile  of  predicted  scalar  dissipation  rate  of  Liew  et 
al.  (1984).  Stretched  flamelet  model  at  x/d  **  20  for  the  values 
of  =  40,100  and  1000. 

Badial  profile  of  predicted  probability  of  combustion  based  on 

stretched  flamelet  model  at  x/d  =  20  for  three  values  of  X  . 

4 

Badial  profile  of  measured  Favre  averaged  mixture  fraction  at  x/d 

=  20  compared  with  predictions  of  model  for  three  values  of  X  . 

4 

Badial  profile  of  measured  mean  density  at  x/d  =  20  compared  with 
predictions  of  model. 


Fig.  21  Badial  profile  of  measured  mean  temperature  at  x/d  “  20  compared 

with  predictions  of  model. 


Fig.  22  Badial  profile  of  measured  Favre  averaged  mixture  fraction  at  x/d 

“  50  compared  with  predictions  of  model. 

Fig.  23  Badial  profile  of  measured  mean  density  at  x/d  -  50  compared  with 

predictions  of  model. 


Fig.  24  Radial  profile  of  measured  mean  temperature  at  x/d  =  50  compared 

with  predictions  of  model. 
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Fig.  14  Meaanred  temperature  ▼».  mixture  fraction  acattergraa  in  the  Ee 

15.000  LHC  flute  at  i/d  -  10,  r/a  -  1.89. 
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Axial  variation  of  aaan  taaparatnra  at  j  -  0  in  tha  Ba  *  1 
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